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Abstract 

The gravitational mass-shift effect is investigated in the framework of the standard model with the energy 

^N . cutoff regularization both for stationary and non-stationary backgrounds at the one-loop level. The problem 

^^ ■ concerning singularity of the effective potential for the Higgs field on the horizon of a black hole, which was 

^ _ reported earlier, is resolved. The equations characterizing the properties of the vacuum state are derived 

, ^ , and solved in a certain approximation for the Schwarzschild black hole. The gravitational mass-shift effect 
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X: 



is completely described in this case. The masses of the massive particles in the standard model are shown 
to depend on the value of the Higgs boson mass in the flat spacetime. If the Higgs boson mass in the flat 



fi \ spacetime is less than 263.6 GeV, then the mass of any massive particle approaching a gravitating object 

I , grows. If the Higgs boson mass in the flat spacetime is greater than or equal to 278.2 GeV, the masses of 

biX) all the massive particles decrease in a strong gravitational field. The Higgs boson masses lying in between 

these two values prove to lead to instability, at least at the one-loop level, and so they are excluded. It turns 

^ ■ out that the vacuum possesses the same properties as an ultrarelativistic fiuid in a certain approximation. 



The expression for the pressure, the entropy and enthalpy densities of this fluid are obtained. The sound 
speed in this fluid is also derived. 
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I. INTRODUCTION 

Recently, many attempts have been undertaken to verify the prediction of general relativity 
in its classical formulation. One of these important predictions is the gravitational redshift law. 
Despite that it was directly tested with a high accuracy [|||, several programs are now planned to 
check the redshift law in the weak field limit to find possible smaller deviation Q. Small deviations 
from the standard redshift law arise in the models that do not preserve the local Lorentz-invariance. 
The latter is believed to be violated on the Planck scale. In this paper we are about to investigate 
a somewhat different way of possible violations of the standard redshift law. It results from the 
gravitational mass-shift effect by means of the Higgs mechanism [^. 

The standard derivation of the redshift law in the course of general relativity heavily relies on 
the assumption that the emission spectrum does not depend on the external gravitational field 
in the reference frame associated with the emitter. However, if the masses of particles (say, the 
electron mass) change with gravitational field, the spectrum will also change even in this system 
of coordinates. The masses of all massive particles of the standard model are generated via the 
Higgs mechanism and determined by a nonzero vacuum expectation value of the Higgs field. This 
expectation value provides the minimum to the effective potential, while the shape of this potential 
turns out to depend on the background gravitational field [^. This leads to the gravitational mass- 
shift effect and, as a consequence, to deviations from the standard redshift law. This is a purely 
quantum effect and it stems from the dependence of the zero-point energy of quantum fields on the 
background gravitational field. 

One of the crucial steps in evaluation of the zero-point energy is a correct definition of the 
Hamiltonian of the whole system in the Hilbert space. The formal expression for the Hamiltonian 
diverges and needs an appropriate regularization procedure. Another equivalent formulation of 
this problem is that we should define the normal ordering for the composite operators y, g]. For 
a stationary gravitational field, we have a distinguished set of the creation-annihilation operators 
which are associated with the stationary mode functions. The latter are defined as eigenfunctions 
of the Lie derivative along the Killing vector related to stationarity (for details, see, e.g., [y|). This 
allows us to define the normal ordering and prescribe a rigorous meaning to composite operators 
such as the Hamiltonian. As soon as the well-defined Hamiltonian is given, the vacuum state is also 
defined as the state with minimum energy (in black hole physics this vacuum is called the Boulware 
vacuum |^). This approach is equivalent to the so-called physical regularization Q, which, loosely 
speaking, consists in modification of the particles dispersion laws above a certain cutoff energy in 



such a way that the loop diagrams become convergent. Both regularization procedures require a 
KilHng vector. It is the square of this Kilhng vector that enters the effective potential of the Higgs 
field, changing the shape of the potential in the presence of external gravitational field. For the 
flat spacetime, this Killing vector "disappears" from the effective action (as its square is a mere 
constant) and, in the one-loop approximation, the effective potential takes a familiar form of the 
Coleman- Weinberg potential [0|. So, after renormalization, there is no Lorentz-invariance violation 
for the flat spacetime according to this approach. 

In the non-stationary case, we have to introduce some vector field ^'^ that allows us to define 
the physical regularization (the notion of a dispersion law) and the normal ordering just as in the 
stationary case. This vector field must coincide with the Killing vector in the case of a stationary 
background. It turns out that this vector field is uniquely defined by the requirement of covariant 
divergenceless of the matter energy-momentum tensor, provided that the system starts its evolution 
from some stationary state. The equations describing the evolution of ^'^ have a hydrodynamic 
form. In particular, there exists a conserved charge, which can be interpreted as the entropy of the 
whole system. As we shall see, in a reasonable approximation, these equations can be cast into the 
form of the Euler equations for an ultrarelativistic fluid with the equations of state determined by 
the effective potential of the Higgs field. This will allow us to find stationary solutions to these 
equations for the Schwarzschild black hole. There are two such solutions. The first one is the trivial 
solution coinciding with the Killing vector. It corresponds to the case of a stable star in a certain 
approximation, of course. The second solution is nontrivial and describes an "accretion" of the 
vacuum onto the black hole. In this case the effective potential turns out to be non-singular on the 
horizon of the black hole. The masses of massive particles of the standard model acquire a finite 
shift and do not vanish there. In general, the gravitational mass-shift effect is greater for stable 
stars than for black holes of the same mass at the same distance from the gravitating object. 

As long as the vector field ^^ enters the effective potential of the Higgs field, it spoils the so- 
called local position invariance (see, e.g., |p!o| ). Hence, the approach we are going to study can 
be somewhat considered among the Lorentz-invariance violating models. Let us emphasize the 



distinctions with the standard methods of introducing the Lorentz-invariance violation |11-15|. 
The vector field discussed above must always be a background field, i.e., it is not a quantum field 
or its average. It characterizes the regularization procedure and necessarily appears in the course of 
defining the composite operators in the Hilbert space. Whatever new particles are introduced into 
the model, this vector field has to be included to define the physical regularization and accomplish 
the theory. The equations of motion of ^^ do not follow from the action principle and arise as 



a self-consistency condition of the model. This approach is minimal in the following sense: As 
we have already mentioned, the Lorentz-invariance is not violated for the flat spacetime after 
renormalization; There is no additional degrees of freedom or structures in the model. In the 
stationary case, the vector field S^^ is a time-like Killing vector and determined by the metric. In 
the non-stationary case, this vector field is uniquely defined by the self-consistency condition. This 
minimality will allow us to describe the gravitational mass-shift effect and make certain predictions 
rather than just to fit various parameters of the model to the experimental data. 

The paper is organized as follows. In Sec. ||, we introduce a notation and recall the basic 
features of the standard model that are necessary for deriving the effective action. As we restrict 
ourself to the one-loop approximation, we need only the spectrum of masses of the real and fictitious 
particles of the standard model. It is obtained in this section for the Feynman gauge. 

In Sec. p^ , we derive the one- loop effective potential of the Higgs field on a stationary background 
in a certain approximation. It has the same form as the effective potential of the Higgs field on 
the Schwarzschild background for the 0{N)-(/)'^ model [^. The different regularization prescriptions 
and the dependence of the final result on them are also discussed in this section. The form of 
the effective potential obtained is completely fixed by imposing three normalization conditions. 
Further, the gravitational mass-shift effect on a stationary background is considered for different 
values of the Higgs boson mass. As expected, there are two scenarios of the mass behavior at small 
the Killing vector squared. The first one corresponds to the infinitely growing masses of massive 
particles at small ^^. The second scenario describes the decreasing masses. In the latter case, 
the symmetry of the standard model is restored at sufficiently small the Killing vector squared. 
Realization of these cases depends on the sign of the undetermined constant resulting from the 
regularization procedure. The sign of this constant depends, in its turn, on the value of the Higgs 
boson mass in the fiat spacetime. This sign changes at the "critical" value of the Higgs boson mass 
which is approximately equal to 278.2 GeV. If the Higgs boson mass is less than this critical value, 
the first scenario is realized. In the opposite case and, in particular, at the critical value of the 
Higgs boson mass, the second scenario takes place. 



Section IV is devoted to the case of a non-stationary background. It begins with derivation of 
the general equations describing an evolution of the vector field ^^ discussed above. Then, in a 
certain approximation, these equations are reduced to the Euler equations for an ultrarelativistic 
fluid. The temperature of this fluid is proportional to {^'^)~^''^, i.e., to the Tolman temperature pq| 
on a stationary background. Therefore gravity effectively heats such fluid. Further, the speed of 
sound for this fluid is introduced. It characterizes the velocity of propagation of small disturbances 



of the vector field (^^. In the weak field limit, this sound speed is independent of particular details 
of the standard model and its square equals the Newtonian potential with reversed sign (the first 
cosmic velocity). The limit of the sound speed at small .^-^ is also universal and equals l/\/3, as it 
should be for a fluid at the very high temperature, although this situation seems not to be realized. 
In a general case, the entropy of this fluid, or the pressure, or the speed of sound cannot be found 
analytically for an arbitrary ^^ in terms of elementary functions. However, if the Higgs boson mass 
is equal (or close) to the critical value, this can be done with a high accuracy for those ^^ that are 
realized in practice. The formula for the sound speed proves to be independent of the details of the 
standard model in this case too. In the interval of the Higgs boson masses from 263.3 GeV to the 
critical value 278.2 GeV, the square of the speed of sound becomes negative for certain ^^ and the 
system is hydrodynamically unstable there. We exclude this region of the Higgs boson masses from 



our subsequent considerations. Section IV is concluded by consideration of the stationary solutions 
to the equations of motion of the vector field ^^ on the Schwarzschild background. Noteworthily, 
the problem of a spherically symmetric accretion admits of an analytic solution for the Higgs boson 
masses close to the critical value. 

We use the metric tensor g^i, with the signature —2 and the system of units in which h = c = 1. 
The Greek indices are raised and lowered by this metric. The gravitational field is assumed to be 
fixed and non-fluctuating. An inclusion of its quantum fluctuations at the one-loop level does not 
influence the results. 



II. ONE-LOOP CORRECTION TO THE EFFECTIVE ACTION 

Let us briefly recall some basic features of the standard model that will be relevant for our 
consideration (see for details [p^7|-|T9[). The generally covariant action of the standard model can be 



divided into several parts 

'S'sM ~ / '^ X^/\g\CsM, -^SM = -^gauge fields + -^leptons + -^quarks + -^ Yukawa + -^Higgsj (1) 

where the Lagrangian densities for the gauge flelds and leptons are 



''-gauge fields — ^JfiiyJ j^^fiu^a a^ fiv'^ a i 
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■-leptons 



i=l 



lliid + |iV, - ^B)ll + lUid - g'B)l^n 



(3) 



Here l]^ are the left-handed SU{2) doublets, /^ are the right-handed SU{2) singlets, and gs, g and g' 
are the coupling constants. The covariant derivative acting on the spinors is defined in the standard 
way 



1 



d = r{d, + -cOf.abh^l'']), 7 



I'e^a, 



(4) 



with Ca and uj^ab being the tetrad and spin connection, respectively. The quark sector reads 
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+ u'nitd + |A-^„ + ^i?)u^ + d'kiid + |A-yl„ - |S)4 , , 



(5) 



d" 






U 



^CKM" ' " -iu,c,t), d' = {d,s,b), 
where f/cKM is the Cabibbo-Kobayashi-Maskawa matrix. The masses of the real massive particles 
of the standard model are generated by the terms 



>Cyi 



"ukawa 



y^W' 
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bn'n + ft 
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<t>cd'n + niWn + h.c. 



(6) 



where /j are the Yukawa couplings and 

■L-Higgs 



k^^M + l^^a + ^i?^; 



]' - ^\^\' - j\<P\\ 



(7) 



21V p ■ 2 '^ " 2 

where (j) is the SU{2) doublet, //^ is a negative constant, and A is the Higgs self-interaction coupling 
constant. 

We consider the minimal standard model, that is, the neutrinos entering into the left-handed 
SU{2) doublets /^ are assumed to be massless and the possible non-minimal coupling ^i?|(/)p of 
the Higgs field with gravity is set to zero. As we shall see, the particles with small masses do 
not considerably change the effective potential and, consequently, do not affect the gravitational 
mass-shift effect. As for the non-minimal coupling term, it can be combined with the term /U^|(/)p 
on the Einstein spaces pO|| . So, as long as we neglect the back-reaction of the matter on the 
metric, this term can be omitted. It should be also noted that if we allow for the gravitational 
field g^j^i, to fluctuate and integrate out the gravitons, the effective Higgs mass is changed by a 



quantity that depends on ^ (see, for example, |21, 22|). However, on the electroweak scale, this 



mass-shift is considerable only for the enormous values of the non-minimal coupling ,^, where the 
very applicability of the perturbation theory becomes questionable. Hence, we disregard this term. 
We shall study the effective action of the Higgs field under the assumption that the vacuum 
expectation values of other fields are zero. To this end, we parameterize the Higgs doublet as 





V + X 



+ iCra 



C2 + Ki 

V + X-iCs 



(8) 



where (^a are the Goldstone bosons, r/ is a vacuum expectation value of the Higgs field, and x 
describes the fiuctuations of the Higgs field. Then the in-out effective action takes the form 



J^g^i^,v] 



/)<J)/)cDPg«'S'tot[9M'-.^:*.C,^] 



IPI 



(out|</>|in) 
(out I in) 



(9) 



'-'tot ^ '-'SM + '-'gauge fixing + '-'ghosts; 

where $ is the full set of fields including the Goldstone bosons and x- while c and P are the 
Faddeev-Popov ghosts. The notation for the path-integral above means that only the one-particle 
irreducible diagrams are taken into account. The states |in) and |out) denote in- and out-vacua. 
They differ only by a phase for stationary backgrounds. In the Feynman gauge we shall use, the 
gauge fixing and ghosts actions are given by 



1 



1 



''-gauge fixing — 9 V^ ^/J ridV^, 



a\2 



Cai^"^ + \9^rf)P'' + c(V2 + jg''^rf)P + c^V^P" + vertices. 



(10) 



■'-ghosts 

where V^ is a covariant derivative and "vertices" denotes the terms of a higher power in the ghost 
fields. These terms do not contribute to the one-loop effective action. 

Now it is not difficult to deduce the spectrum of real and fictitious particles of the standard 
model in the Feynman gauge. It is presented in Table |. This is the only information about the 
standard model that we need for the one- loop calculations. Notice that all the masses of massive 
particles are proportional to r/ except for the Higgs and Goldstone bosons. Furthermore, small 
deviations of ry result in the linear response of all the masses. 

The Killing vector ^^ characterizing stationarity of the system can be straighten (^^ = (1, 0, 0, 0)) 
by an appropriate coordinate changeQ. Then the formal expression for the one-loop contribution of 
one bosonic mode to the effective action on a stationary background can be cast into the form 



16 



k 



(11) 



^ We use the standard definition for stationary and static spacetimes (see, e.g., [E4|). Tlie spacetime is called 
stationary if it possesses a Killing vector ^^ such that ^^ — g^.v^^C' > 0. The stationary spacetime is said to be 
(locally) static if V[^(^^]/,f^) = 0. 



Name 


Mass squared 


# 


Higgs 


SAry^ + ;,2 


1 X 1 


Goldstone Z 


TO^ +m,| 


1 X 1 


Goldstone W 


m^ +m^ 


2 X 1 


Z 


(g2+g'2)^2/4 


1 X 4 


w 


5'r;V4 


2x4 


Ghosts to SU{2) 


m^ 


-6 X 1 


Photons 





1 X 4 


Ghosts toU{l) 


m|-m2^, 


-2 X 1 


Gluons 





8x4 


Ghosts to SUiS) 





-16 X 1 



Name 


Mass 


# 


Ui 

d. 


1' 


9x4 
9x4 


e,fj,,T 


3x4 


H 





3x2 



Quantity 


GeV 


^70 


247 


mm 


129 


mzQ 


91.2 


mwQ 


80.4 


mta 


171.2 


rribo 


4.2 


rrico 


1.27 


rriro 


1.777 



TABLE I. The spectrum of particles of the standard model in the Feynman gauge. The bosons and ghosts to 
them are presented in the left table. The characteristics of the fermionic particles are given in the center table. 
For brevity, we introduced the notation ffi^ :— \rf + /i^. The last column in these tables contains the number of 
degrees of freedom for each particle: the number of particles A'^ of a given type multiplied by the number of their 
polarizations in the Feynman gauge. The ghost fields have a negative number of degrees of freedom by definition. 
The right table collects the vacuum expectation value of the Higgs field and the masses of particles that are greater 
than 1 GeV. The experimental data were taken from l23] and correspond to the renormalization scale 2 GeV. The 
recommended value of the Higgs boson mass is also given. 



where T is the time interval and E^^ is the energy of the mode k corresponding to an antiparticle. 
This contribution is proportional to the zero-point energy and results from the normal ordering of 
the creation-annihilation operators. Another equivalent form for the one- loop contribution reads as 



ri;^ = -TSp r ^^(-G-^(-po)) = -TImSp /~ |^ln(G-( 



-Po)), 



(12) 



where G{po) is the Fourier transform in x^ of Green's function corresponding to the field considered 
on a curved background. The Heaviside step function 6{L — A) of the Hermitian operator L defines 
the spectral decomposition of unity associated with this operator. As for fermions, they give the 
same contributions but with an opposite sign in Eq. (0). The complete one-loop contribution to 
the effective action of the standard model on a stationary background is given by the sum of the 



contributions (11) over all the species of particles multiplied by their number A^ presented in Table 



III. EFFECTIVE POTENTIAL 



Now we are in position to calculate the one-loop effective action. To this end, we have to 



prescribe a rigorous meaning to the formal expressions (O) or (12) and regularize them. One can 



distinguish two ways, at least, how to do this: we can employ the commonly used regularization 



methods based on the notion of a Fock proper-time |25-27|, or introduce the energy cutoff - the so- 
called physical regularization (see, e.g., [P|). The first approach uses the heat kernel representation 
for the complex power of the Feynman propagator 

GF = -f^^^i^s)'-'e-^-\ A€C. (13) 

As usual, the Feynman propagator Gp is specified by the mass shift m? — t- m? — ie. The one- 
loop correction can be expressed in terms of the trace of the operator G^ (see, e.g., p8|-p0|). The 
regularization is achieved by the analytic continuation in A (the zeta- function regularization), or 
in the parameter d ~ the dimension of the spacetime - entering the heat kernel (the dimensional 
regularization), or by a mere cutoff on the upper integration limit. All these regularizations, which 
we shall collectively call the proper-time regularizations, give the same answer up to redefinition of 
the arising infinite constants. This is the general property of regularization schemes depending on 
one parameter. They are all equivalent in the aforementioned sense provided the multi-dimensional 
divergent integral (the trace, in our case) is reduced to the one-dimensional divergent integral over 
the same variable (the proper-time s, in our case). The regularization of one-dimensional integrals is 
uniquely defined up to the freedom mentioned above pl| . As regards the multi-loop contributions, 
in the Feynman parameterization, they are all reduced to the one-dimensional integral with respect 
to the sum of the proper-times. 

The second regularization prescription can be conveniently realized by introducing the Fermi- 
Dirac thermal cutoff. Namely, consider the one-loop fi-potential for fermionic fields at zero chemical 
potential (in the system of coordinates where ^^ = (1,0,0,0)) 

o «-i, yfR^ Q r^ e{G~H-po)) [^ dpo ln{-G-\-po)) 

n = P, lnZ(/3o) = Spy^ dpo ^^^^^^^ =^"^^Pyo ~ e^opo + i ' (1^) 

with /3o being a reciprocal temperature defining the cutoff. Then the average energy takes a familiar 
form and reproduces the zero-point energy in the high-temperature limit 

E = -dsJPo^) = y ^ >y^. (15) 

This representation allows us to employ the known high-temperature expansions for the partition 
function of fermions to obtain the one- loop effective action. Other regularization schemes intro- 



ducing the energy cutoff into the model will lead to the same answer for the effective action up to 
redefinition of the infinite constants. 

In this paper, we shall investigate the second method of regularization - the energy cutoff - for 
the following reasons (see also ||8|). The energy cutoff^ has a clear physical interpretation in terms of 
particles: the dispersion laws of particles just deviate from the relativistic dispersion law at very high 
energies. This method is also natural from the quantum mechanical viewpoint as it is formulated in 
the quantum mechanical terms. It does not lead to states with a negative norm. This regularization 
procedure is minimal, i.e., we need not to introduce new entities into the model (as, for example, 
the proper-time) when the notions of energy and dispersion law are defined. However, this method 
possesses the well-known drawbacks: it violates the Lorentz- and gauge symmetries and requires 
a fine-tuning of the coefficients at the counterterms to restore these symmetries in the effective 
action (see, e.g., pi p^); it is awkward in comparison, for example, with the standard proper-time 
dimensional regularization. As opposed to the energy cutoff, the regularization procedures using 
the proper-time representation for the Feynman propagator preserve almost all the symmetries of 
the model since the proper-time is a Lorentz- and gauge invariant quantity. However, the concept 
of proper-time is purely classical and it cannot be defined in quantum mechanical terms without 
enlarging the number of degrees of freedom (introducing, for instance, one extra dimension for it). 

On a flat background, these two approaches to regularization are equivalent (see, e.g., ||8|, |32| - 
|35| ) in the sense that we can adjust the coefficients at counterterms so as to make the effective 
actions equal. However, as we shall see, these two schemes (^) and ([iq ) give inequivalent results 
on a curved background. The divergent and finite parts of the effective actions will differ. At 
that, as for the finite part, this inequivalence is nonpolymomial in fields and cannot be removed by 
counterterms or field redefinitions. 

The high-temperature expansions of one-loop contributions to the fi-potential on an arbitrary 
curved background are rather huge (see, e.g., p6|). That is why we restrict ourself to finding 
only the effective potential of the Higgs field. We neglect all the derivatives of the Higgs field. 
Furthermore, we also discard all the terms proportional to the derivatives of the metric field. The 
former approximation is adequate when the Higgs field changes slowly on the scales characterizing 
the metric field, while the latter approximation is appropriate when the background metric is 
approximately Ricci-flat. In this case, most of the relevant terms of the high-temperature expansion 



■^ Strictly speaking, we need to cut out a finite dimensional subspace from the one-particle Hilbert space. Only in 
this case do we remove all the divergencies. To this aim, we can confine the particle into a sufficiently large "box" 
making its spectrum discrete and then single out the subspace in the one-particle Hilbert space. This subspace 
spans the mode functions with the energies not higher than the energy cutoff. 
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proportional to the derivatives of the metric vanish. The only possible nonvanishing combinations 
containing derivatives and having appropriate dimensions are (the total derivatives are also omitted, 



cf. Pi, Eqs. (4.27), (4.28)) 



V.lnVi^V^lnV^, 



R. 



/lupcr 



R^'"^^, (V^lnx/i^Vnnv^)2, Vu(e.i/e')V[^(rVe'), (16) 



where ^^ = Qixv^'^i^ ■ The last term in (^) vanishes for static metrics. We discard these terms and 
shall obtain thereby the leading contribution to the derivative expansion of the effective action. As 
we shall see, this approximation is quite reasonable for the metrics of stars and macroscopic black 
holes. Every derivative of the metric or Higgs field diminishes, effectively, the contribution of the 
term to the effective action by the factor (mr)^^, where the massive parameter m is of the order 
of 1 — 100 GeV (electroweak scale) and r is a distance from the gravitating body. 

Bearing these assumptions in mind, it is easy to calculate the O-potential (14). The simplest 
way to do this is to use the representation for the trace of an operator in terms of its symbol (see, 
e.g., |38-fll[). If we neglect all the derivatives of the fields, then the symbol of the operator entering 



(14) is given by 



for a scalar field. Hence, the partition function takes the form 



(17) 



In Z(/3o) f« /3o / '^^0/7^ 



^ r ^ f dxdp e ( (5°° - gijo^'g'^^ )pI + g'h^P3 - m^) 
f3o dpo ' ^ 

/3o / dpo 



(27r)3 e*Po + 1 

dxdp 6 {pI + gooig'^PiPj - w?)) 



(18) 



(27r)3 



(2vr)3 V500 



e/3opo + 1 
In ( 1 + e"^"^ooVp^+™^') 



where we have used the relations p2l 43 1 



^^ = 6i , 500 det gi, = g, 5°° - 5^,5°*5°^' = (500)"^ 



gikg 



Thus, the contribution of one bosonic mode to the effective potential has the form 



.(1) 
lb 



T / dxyyids 



dp 

(2^ 



In (1 + e-/^Vp'+™' 



(19) 



(20) 



/3o-^0 



where /3 := /3o\/^ is the Tolman reciprocal temperature. The contribution of the spin degrees of 
freedom just multiplies ( pO| ) by the corresponding factor. We have derived the well-known result 
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1 16, Eg! that, in the leading quasiclassical order, the influence of gravity on thermal contributions 



reduces to the substitution of the temperature by the Tolman temperature. Now we can use the 



high-temperature expansion of the Fermi-Dirac il-potential on a flat background (see, e.g., |18 



r« = T / dxVl5l 



+ 777-0 In ^^0_^ + 27+- 



, (21) 

/3o->0 



240/3^^4 48/3^^2 647r2 V 7r2 ' ' ' 2 

where 7 is the Euler constant. The same answer would be obtained if we used the point-splitting 
regularization procedure (|P, Eq. (247)). The one-loop correction (^) to the effective potential 
is generally covariant and represents the leading term in the derivative expansion of the effective 
action. The Killing vector entering the effective potential is not an external structure on the 
spacetime. It is determined by the metric through the Killing equations. Formula (^) on a 
spherically symmetric background can be also derived using the explicit quasiclassical representation 
of the mode functions |^. The total one-loop contribution to the effective potential is obtained by 
summing over all particles. 

If we employed the proper-time regularization we would obtain the same result (21), but with 



^2 replaced by one. This fact is quite expectable as the rate of the proper-time of a particle moving 
along the Killing vector ^'^ is proportional to y^- 

ds = y^dt, (22) 

where t is the variable dual to the energy (time). Loosely speaking, when we use the energy 
cutoff, the regularization parameter becomes blue-shifted (in comparison with the proper-time 
regularization) by the gravitational fleld. This observation suggests a simple way how to modify 
the results obtained by means of the proper-time regularization to the energy cutoff regularization. 
An accurate inspection of this observation will be given elsewhere. 

So, as we see, there are three type of divergencies in the effective potential. They have to be 
canceled by adding appropriate counterterms to the initial action of the standard model. This, 
in essence, reduces to replacement of the infinite constants (3^ , (3q , and ln/3o by some finite 
quantities; and all other quantities entering the effective potential are replaced by their renormalized 
values. After that, we should impose certain normalization conditions in order to fix five unknown 
constants A, yU^, A, ^O) and A2, where A is related to the cosmological constant, while ^0 ^ind A2 
are the dimensional constants at ^"^ and S,'"^, respectively. 

The first normalization condition we choose is that the effective potential reduces to the well- 
known Coleman- Weinberg potential 0] for the fiat spacetime, i.e., at ^^ = 1^ and the masses of 
particles take their experimental values there. It is convenient to take the values of masses on the 
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scale 2 GeV (see Table ||) where the standard model is in perturbative regime and the one-loop 
approximation is reasonable. This normalization condition implies 



V = A + ^r]^ + -rj 



> rrif In -^ > m? In — 2=— 



2 4 647r2 , ^ „.^Q ^ „.^Q 






(23) 



where the indices / and h indicate statistics and ttiq's denote the values of the masses of the particles 
at ^2 = 1. The coupling constants [i} and A are determined by the equations 



/^^ 



-A^o + :7^{X]"/"^/o-Z]"''"^w)> "T-Ho = 3Ar/g + ^^ (24) 

where a's are defined by brr? =: adrf'. This system of equations can easily be solved exactly, but 
its solution is rather awkward. For a rough estimate, the term in parenthesis in the right-hand side 
of the first equation can be omitted. 

As for the second normalization condition, we demand that the energy-momentum tensor follow- 
ing from (^) vanishes at .^^ = 1 (henceforward, we neglect the cosmological constant). It follows 
from this condition that 



W/^2 2^1 _n dVie^r,^) 



\e=i "' de 



0, (25) 



where r] is taken at the minimum of the effective potential V. These equalities fix the constant A 
and relate the constants Aq and A2'- 

^ = -y ^0 - 4%'> A (E - E) = "I^^ (E "^/o - E "^'oj - y (E "^/o - E ""boj ■ 

(26) 
Thus, we ascertain in our particular case that the energy cutoff and the proper-time regularizations 
give equivalent results for the effective action on a flat background (^^ = 1, in our case) and are 
inequivalent in a curved spacetime. 

The only undetermined constant can be fixed imposing the third normalization condition. 
Namely, the deviation of the vacuum expectation value ry due to change of ^^ leads to the de- 
viation of masses of all the massive particles and, as a consequence, to the deviation from the 
standard redshift law. Recall that in general relativity the redshift law is usually derived under the 
assumption that the emission spectrum does not depend on gravity in the reference frame associ- 
ated with the emitter. However, if the masses of massive particles (for example, the electron mass) 
change with the gravitational field, the spectrum will also change even in this system of coordinates. 
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FIG. 1. The dependence of the constant A2 on the Higgs boson mass according to the normahzation condition 
(ESI). Thickness of the hne depicts the error related to the unknown constant a. Notice that A2 is almost uniquely 
determined and can be approximated by its value at a = 0. The inset represents the dependence of the Higgs 
self-interaction coupling constant on the Higgs boson mass. 



It should be noted that the null redshift experiment |10| cannot apparently detect the mass-shift 
effect, at least in the linear order in gravitational perturbations. The additional corrections due to 
mass-shift to the emission spectra are the same as the changes of the linear sizes of the resonator 
in superconducting-cavity stabilized oscillator clocks. The sizes of the resonator vary in the same 
way as the Bohr radius of the hydrogen atom. So, the direct experiments are necessary. The direct 
tests of general relativity impose stringent constraints on the deviation from the standard redshift 
law in the weak field limit ||l], Q: 



UJ 



{l + a)LpN, 



lal <2 X 10"^, 



(27) 



where (/?7v is the Newtonian potential and uj is the frequency of radiation. It is not difficult to 
brought this condition into the form 



5rl_ 



adf, 



e=i 



(28) 



where rj is taken at the minimum of the effective potential as before and 5S,'^ ~ 2^9 a^- This condition 
fixes the last undetermined constant in the effective potential (^). Making a negligible error, we 
can set a = (see Fig. Il|). Then 



dV{^'^,r]^ 



Q^2Qj^2 







A2 



1 Ef"/ 



m 



/o 



Efe «fe"^bo 



327r2 



(29) 



«=1 -J^"- E/"/-E6"b 

As a result, we can completely describe the gravitational mass-shift effect in the standard model 
on a stationary background. 

The standard model is not an accomplished theory. The Higgs boson has not yet been found and 
its mass is a free parameter of the model. Therefore several possible scenarios are presented in Fig. 
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^. They are the same as for 0{N)-(j)'^ -model considered in Q, but "inverted", because the fermionic 
degrees of freedom dominate in the standard model. There is a critical value of the Higgs boson 
mass wi^o ~ 278.2 GeV where the constant A2 vanishes. When the Higgs boson mass is equal to 
or greater than m'^Q, the masses of all the massive particles decrease with the gravitational field 
and, eventually, go to zero at sufficiently small ^^. The broken symmetry of the standard model is 
restored. The system passes through the phase transition of the second order for mno > '^i/o- ^^ 
case ruHo = ''tT'hoi it passes through the phase transition of the first order and then of the second 
order (see Fig. g). When the expectation value rj decreases, the masses squared of the Goldstone and 
Higgs bosons become negative. This leads to an appearance of the imaginary part in the effective 
potential and says that the vacuum state gets unstable. The values of the logarithms entering the 
effective potential taken on their cuts are uniquely specified by the prescription m^ — )• m? — it. 
Strictly speaking, in this case we have to use the in-in formalism for systems with a non-stationary 
vacuum state. However, if the imaginary correction to the expectation value is rather small. 



ImV;',(e2,r/2 



Re(r?2y'^^,(^2,^2)) 



< 1, (30) 



we can still use the in-out effective action for the approximate evaluation of averages. The rate of 
decay of the vacuum state P7| , i.e., the probability, per unit time per unit volume, that a pair is 
created, is approximately given by the imaginary part of the effective Lagrangian doubled. In our 
case, it is nonncgative, as it should be, and is written as 

-2Imy = ^ [m\je{-ml) + 2m%we{-ml^) + m%z0[-mlz)] , (31) 

where niGW a-nd mcz are the masses of the corresponding Goldstone bosons. The imaginary part 
arises when 

W%'<(l + 2m2^oMD"'- (32) 



The numerical analysis shows that condition ( |30|) is fulfilled for any ^^ at which the phase transition 
occurs and, of course, for larger ^^ (gee Fig. |2|). 

In the opposite case, i.e., when rriHQ < m'^^, and, in particular, for the recommended value of 
the Higgs boson mass mj^Q = 129 GeV, the masses of massive particles grow, when ^^ decreases, 
and tend to infinity at ^^ — )• 0. At small ^■^, the asymptotics of the expectation value of the Higgs 
field is given by 

7?Vr?o' ~ c/f, (33) 
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FIG. 2. Left top panel: The dependence of the vacuum expectation value of the Higgs field on ^^ at the recommended 
value of the Higgs boson mass. The left inset represents the same dependence, but at smaller ^^. The right inset 
depicts the relative error 1 — r\aslr\ of the asymptotics (631). Right top panel: The dependence of the vacuum 
expectation value of the Higgs field on Ig^^ at mno = 279 GeV. The left inset represents the dependence of the 
imaginary part of the Lagrangian doubled. The right inset provides the dependence of the quantity (bol) on Ig^ . 
This quantity characterizes the applicability of the in-out formalism in the case at issue. Left bottom panel: The 
same as for the right top panel, but for the critical value of the Higgs boson mass m^g ~ 278.2 GeV. Notice that 
the masses do not almost change in this case provided ^^ is not extremely small. Right bottom panel: The shape of 
the real part of the effective potential is shown in a vicinity of the phase transition for the critical value of the Higgs 
boson mass (the value of the real part of the effective potential at 77 = is subtracted). The solid line corresponds 
to Ig^^ = —16.785, the dotted line is for Ig^^ — —16.75, and the dashed line represents the case Ig^^ = —16.8. It is 
clearly seen from these plots that the first phase transition is of the first order. 

where the dimensionless constant c is found from the equation 






IA2 



+ A 






al 



1 + In c^ 

327r2 



1 



167r2 



a]j In ^^^ + 2a%^r In 



aawVo 



m. 



GZO 



aczrio 



c = 0. (34) 



It has the solution c ~ 9.75 x 10"^ at the recommended value of the Higgs boson mass. In Fig. ^, 
this asymptotics is presented in comparison with the exact expectation value. 

Some comments on these scenarios are in order. At small ^2, the quantum correction becomes 
large and the gradients of the fields ^^ and 77 increase too. Therefore we need to take into account 
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the higher-loop corrections, the gradient terms, and back-reaction when ^"^ tends to zero. Moreover, 
the Higgs field self-interaction coupling is large at large the Higgs boson mass tuhq'- A ~ 0.63 at 
the critical value of the Higgs boson mass. Hence, the critical value given above is just an estimate, 
and the scenarios described for large rriHQ are rather qualitative than quantitative. The lattice 
simulations and higher-loop calculations for the standard model at finite temperature nevertheless 



show that such one-loop results are quite reasonable (see, e.g., |18, pH|). In the next section we 
shall see that very small values of ^^ seem not to realize even for black holes. Also notice that the 
effective potential (|23| ) evidently depends on the spectrum of masses of particles of the standard 
model. It will change if some new heavy particles are discovered. 

IV. NON-STATIONARY CASE 

In the previous sections, we considered the standard model on a stationary background when the 
Killing vector ^'^ is defined. In the long run, the square of this vector entered the effective potential 
of the Higgs field. A natural question arises how to generalize the above results to non-stationary 
backgrounds. It turns out that this generalization is essentially unique under certain reasonable 
assumptions. 

To investigate the dynamics of quantum field on a non-stationary background we have to use the 



in-in formalism and the in-in effective action (see [g^, 46-49|]). This approach doubles the number 



of the quantum and background fields {g^iy, $) — )• {g^^, ^^)i where ^ = (<^, c, P). Further, we must 
prescribe the exact meaning to the operators entering the Heisenberg equations and introduce an 
appropriate regularization. The self-consistent definition of operators in the Fock space needs the 
definition of the normal ordering [^, ||, 38 1. Different prescriptions for the normal ordering lead 



to inequivalent theories on a curved background (see, e.g., |^), while the formal manipulations 
with divergent operators may give rise to incorrect results like the absence of anomalies (whereas 
they are known to exist). On a stationary background, we have a preferred set (or sets) of the 
creation-annihilation operators and, consequently, a preferred prescription for the normal ordering. 
These creation-annihilation operators are associated with the stationary mode functions that are 
eigenfunctions of the Lie derivative with respect to ^^. In the case of a degenerate one-particle 
energy spectrum, there is a freedom to combine the mode functions corresponding to the same 
energy, but this freedom does not affect the observables. It just results in a unitary transform in 
the finite dimensional subspace of the one-particle Hilbert space provided the particle is confined 
into a large "box". Generalizing this construction to non-stationary backgrounds, we introduce the 
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background vector field ^^ as well and require that this vector coincides with the Killing vector in 
the case of a stationary background. As before, it is normalized on unity for the flat spacetime: 
^^ = 1. This vector field allows us to define the normal ordering and the physical regularization. 
The regularized in-in effective action becomes the functional of the form 

r[5+.,e+,^+;V,e;,'^-]- (35) 

It is generally covariant with respect to the first and second arguments, separately. This implies 

<5r 



v-v « £^r^ + V'i.T^, v^r^" := -2 ^ 



l5|r» := 



6^-^ 



(36) 



Here the approximate equality means that the equations of motion for the fields <I> are taken into 
account and we identify the "plus" and "minus" fields upon variation. We see that covariant diver- 
genceless of the energy-momentum tensor (the average of the energy-momentum tensor operator) 
is violated by the terms depending on the vector field ^'*. In the stationary case, these terms disap- 
pear since ^^ is the Killing vector (see Eq. (|36|)), and the energy-momentum tensor is covariantly 
nondivergent. In the non-stationary case, we just impose this condition on the energy-momentum 
tensor and fix thereby the vector field ^^. Recall that the divergenceless of the energy-momentum 
tensor follows from the Einstein equations and so the latter requirement is a mere self-consistency 
condition (a Ward identity). The equations of motion for the vector field ^^ can be cast into the 
hydrodynamic form 

V^^iew) = 0, e''V[^(r,]/ti;) = C^{Tjw) = 0, (37) 

where w := ^^Tp. These equations should be supplemented by the initial and boundary conditions 
discussed above. It follows from the first equation that there exists a conserved charge in the system. 
Matching Eqs. (p^) with the equations of motion of the relativistic fiuid (see, e.g., ||50|]), it is natural 
to identify this charge with the entropy of the system. The second equation implies that if the 1- 
form uj^^Tfj^ was initially exact, it remains exact along the integral curves of the vector field ^^. 
In this case, the spacetime is foliated by the hypersurfaces associated with this integrable 1-form. 
This provides a preferred definition of the energy of the system on a non-stationary background as 
a flux of the vector field T^'^^y through one of these hypersurfaces. Of course, this energy is not 
generally conserved as long as the background is non-stationary. 

In a certain sense, the vector field ^^ spoils the so-called local position invariance (see, e.g., 
pop. Hence, the approach we are considering can be placed among the Lorentz-invariance violating 
theories (for a review, see p3| , |T^]). Notice, however, a distinction of this approach with the standard 
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Lorentz-invariance violating models like, for example, the minimal Lorentz-invariance violating 
standard model extension ||l3| . The vector field ^'^ we have introduced is always a background 
field, i.e., it does not represent a quantum field or its average. It characterizes the regularization 
procedure and is necessary for a proper definition of the composite field operators in the Hilbert 
space. Whatever new particles are introduced into the theory, this vector field has to be included 
when it comes to define the physical regularization. Therefore its dynamics do not follow from 
the action principle, at least immediately, but arise as a self-consistency condition. Moreover, this 
approach is minimal: i) There is no Lorentz-invariance violation for a flat background; ii) There 
is no new additional structures or flelds in the effective action on a stationary background. The 
vector field ^^ is just the Killing vector of the metric. In the non-stationary case, this vector field is 
uniquely determined by the self-consistency condition (|37| ) with the initial and boundary conditions 
described above; iii) The number of new parameters entering the effective action is minimal, which 
allows us to make predictions rather than just to fit the model to the experimental data. 

As an example, we consider the Schwarzschild black hole. Despite the fact that it is described by 
a static metric except for the small region near the horizon where the matter accretes, it should be 
thought of as an infinitely lasting collapse, i.e., as a non-stationary system. Later, this observation 



will be relevant for the analysis of the solutions to (37). 

Now we should make certain approximations in order to solve the equation of motion for the 
vector field ^^. First, we assume that the system is in a vacuum state, i.e., there are no particles in 
it, or their back-reaction on the vacuum is marginal, and the Hawking particle production [|5l| is 



also negligible for the averages of quantum fields. Then the in- in effective action entering (37) can 
be approximated by the in-out effective action for the vacuum: 

Second, we use the same assumptions which we made in deriving the effective potential in the 
previous section. That is, we replace the effective action entering (|37| ) by the effective action 
constructed from the effective potential (ESI). Then we have 



r. = 2e.^^qy^, y^ = 2e^-^^^l^, r,/w = ue, (39) 



where r] is taken at the minimum of the effective potential. The last relation in ( |39| ) holds in 
a general case provided that the vector field ^'^ appears in the in-in effective action as S,"^. The 
combination S,^/^"^ is nothing but the Tolman reciprocal temperature 1-form when ^'^ is the Killing 
vector. It is closed when the metric is static and, hence, it is exact provided the fundamental group 
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of the spacetime is trivial. So, if the system starts its evolution from the state with a static metric, 
the 1-form w'^T^ will be exact over all the spacetime in this case. 



The equations of motion (p7|) together with (39) are exactly the equations describing a hot 



ultrarelativistic fluid |50| 



/3:=v^, a:=l3w, p := -V, (40) 

being the reciprocal temperature, the entropy density, and the pressure, respectively, and w =: e+p 
being the enthalpy density (for the other hydrodynamic descriptions of a vacuum, see, e.g., ||52| ). As 
usual, the energy density is denoted by e. The temperature scale is chosen such that the temperature 
(C^)~^ is equal to unity for the flat spacetime. As a result, this temperature is dimensionless and 
the entropy density has the dimension of the energy density. These quantities can be also obtained 
from the standard definition of the energy-momentum tensor if we vary the effective action with 



the potential V{^'^,r]'^) with respect to the metric. The normalization conditions (25) imply that 
the entropy, enthalpy, and energy densities, as well as the pressure of the vacuum state vanish for 
the flat spacetime. 

In the stationary case, we have already seen that the physical regularization can be achieved 
by introducing a thermal cutoff as if the system were at the very high temperature characterized 
by the energy cutoff. Therefore it is not surprising that, in the non-stationary case, a vacuum 
acquires the properties of an ultrarelativistic fluid. According to this interpretation, one can think 
of the particle-antiparticle virtual pairs as some particles created at the expense of the energy of a 
thermostat heated up to the cutoff temperature. This energy is given back to the thermostat when 
these particles annihilate. All the charges are conserved during this process and the average energy 
is also conserved in the stationary case. Of course, this is just one of the possible ways of thinking 
of the vacuum and particle-antiparticle virtual pairs likewise the notion of Dirac's "sea". One may 
say that this sea is heated up to the cutoff temperature. 

Inasmuch as the vacuum possesses the same properties as an ultrarelativistic fluid, we can define, 

in particular, the speed of sound in it 

^2 .^ dp f dlnay' 

' ■ de \d\nf3j ' ^ ^ 

i.e., the speed of propagation of small perturbations of the vector field £^^. Bearing in mind the 



normalization condition (25), we can easily find the expression for this speed in the weak field limit 



ci « ^^ « -V^N, (42) 



„2_,1-C^ 
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where we have assumed that 



a 



ie)L^ 



2p"{e),.^, 


-1 

167r2 











is not zero. The last formula is easily derived if one takes into account the normalization condition 
( p9| ) . It is remarkable that the speed of sound is independent of details of the standard model in 
the weak field limit. It is solely determined by the Newtonian potential ipN and equals the so-called 
first cosmic velocity. One may say that the virial theorem is fulfilled for this fluid. In the case when 
cr'{l) = 0, while cr"{l) ^ 0, the sound speed squared in the weak field limit is given by (H2|) divided 
by 2. This occurs at the Higgs boson mass friHo ~ 263.6 GeV, but it is hard to imagine that this 
degenerate case is realized indeed. It is this case when the components of the energy-momentum 
tensor of the vacuum tend to zero as l/r^ and not as 1/r, where r is a distance from the gravitating 
object. 

At small ,^2, the massless contribution to the entropy density dominates. This is not difficult to 
see from the asymptotics (|3^ ) when the masses of particles grow with decreasing ^"^ . In the opposite 
case {A2 > 0), the masses tend to zero and, of course, the massless contribution also dominates. In 
both cases, 



a^ ^' 



^ (Yl "'/o - Yl ""w j + 2^2 (Y, "^/o - Y ""w j 



327r2 , 

/ b ' -^ / 



TT *^'' 



as one would expect for a gas of massless particles. 

As seen from Fig. |2[ the masses of particles at the critical value of the Higgs boson mass {A2 = 0) 
do not almost change down to extremely small values of ^^ , This allows us to find exact expressions 
for the entropy density and the speed of sound for sufficiently large ^2; 

^= (Y'^'fo-Y'^boj 32^2 ' ^' = ^^7^- ('i^) 

Note that this speed does not depend on the parameters of the standard model. Its weak field limit 
coincides with (42), and it turns into (p4|) at ^^ _), Q. The plots of c^ versus ^^ are presented in Fig. ^ 



for different values of the Higgs boson mass. In the interval of the Higgs masses niHo € {rhno, "^hq)^ 
there is a range of the values of ^^ where the sound speed squared is negative and the entropy density 
changes its sign (see Fig. ^. This signalizes that the system becomes hydrodynamically unstable, 
i.e., small fluctuations of the vector fleld ^^ grow exponentially with time. Since this occurs at large 
values of the Higgs boson mass uiho-, such a behavior may be a mere artefact of the approximations 
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FIG. 3. Left panel: The dependence of the sound speed squared on 1 — ^^. The sohd hne represents the sound 
speed at the recommended value of the Higgs boson mass. The dashed hne corresponds to the critical value of the 
Higgs boson mass rrfno- This plot coincides with the plot of the function (ka). The dotted line depicts the sound 
speed squared at the Higgs boson mass mH=o ~ 264.6 GeV belonging to the interval of hydrodynamical instability 
{niHOy'm'Ho). The inset represents the plot of the relation Cso/cs — 1, where Cso is given by formula (WSl) and Cs is 
the sound speed at the recommended value of the Higgs boson mass. Right panel: The dependence of the entropy 
density of the vacuum state on the reciprocal Tolman temperature squared ^^. The solid, dotted, and dashed lines 
correspond to the same cases as on the left panel. The inset depicts the same dependence, but for smaller 5^. We 
see that the entropy density is positive in the case A2 > and is mostly negative otherwise. The entropy density 
changes its sign with ^^ at the Higgs boson masses lying in the "instability" interval (mHO,"T-Ho)- In particular, the 
entropy density is nonpositive at the recommended value of the Higgs boson mass. 

made and could be cured by higher-loop and derivative corrections to the effective action. However, 
this question needs a further investigation. Henceforth, we assume that the Higgs boson mass does 
not fall into this interval and the system considered is hydrodynamically stable. 

Now we consider in detail the solution to the equations of motion (^) for the Schwarzschild black 
hole. Upon the approximations made, we have just to describe a spherically symmetric accretion 
of an ultrarelativistic fluid onto a black hole. This is a well-studied subject (see, e.g., p3|-|56|) and 
we consider only the main steps. In this case, the self-consistency condition (|37| ) is reduced to two 
equations 



1 



— + U 



r^uo{i' 



k, 



(46) 



where Vg is the gravitational radius, k is some constant characterizing the entropy flux, and u := 
C^/vs is a radial component of the 4-velocity. We have also used the normalization condition 
^^ — )• 1 at spatial inflnity. The system of equations (|4q ) has the trivial solution 

-2 



k = 0, u = 0, 



C 



1-^. 



(47) 



It corresponds to the static system studied in the previous section or, physically, to a stable star (in 
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a certain approximation, of course). Stability of the star means here that it does not appreciably 
change its gravitational field on the time scale of the order of L/cg, where L is a characteristic size 
of the star. A black hole represents an infinitely lasting collapse and so the entropy flux and the 
radial velocity u are not zero in this case. We see that the Tolman reciprocal temperature /3 for a 
black hole is greater than the reciprocal temperature for a static star at the same coordinate r. In 
particular, the reciprocal temperature is not equal to zero at r = r^ for a black hole. Therefore the 
gravitational mass-shift effect is greater for a static star rather than for a black hole at the same r. 
The system of equations ( |4q ) has a unique nontrivial {k 7^ 0) nonsingular solution u{r) satisfying 
the boundary condition 



lim u = 0. (48) 

r— s-oo 

It possesses the weak field asymptotics 

u ^ -X , (49) 

provided muQ 7^ muQ. The constant k is uniquely defined by the regularity condition. A regular 
solution u(r^ should pass through the so-called critical point where, roughly speaking, the velocity 
of the accretion flow becomes equal to the speed of sound in it: 

ul = iyM) = -^^. (50) 

The asterisk distinguishes the quantities taken at the critical point. The equations for the critical 
point follow from (H^) . If we differentiate (H^) then 

(^ + 2nV)— + (2ra + T,a')% = 0. (51) 

The quantity, 

,2 



u 

ec 



a + 2nV = a ( 1 - ^ ) , (52) 



changes its sign when r runs from Vg to -|-cxd since 
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o- 1 - 7:5-^ — > a{l- Cg ), all- —t-t: — > a. (53) 

The entropy density has the same sign for any ^ G (0, 1) and the sound speed is always less than 
unity provided the Higgs boson mass for the flat spacetime does not fall into the "instability" interval 
(mHOi'^T'Ho)- '^^^^ requirement that the radial velocity has a finite derivative at the point, where 
the quantity (J5^) vanishes, leads to equations (50) for the critical point. Equations ( |50| ) allow us 



to find r*, u*, and ^J, and, hence, the constant k from the second equation in (|4^). 
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Unfortunately, the critical point cannot be found analytically with the exception of the case 
where the Higgs boson mass niHo is equal (or close) to the critical mass m'j^Q. However, the 
accretion problem can be easily analyzed numerically. The plots of the accretion velocity and the 
Tolman temperature are presented in Fig. |^for the different Higgs boson masses. The critical point 
for the recommended Higgs boson mass is found to be 

C^ w 0.90, u* w -0.18, r^/rg ^ 7.73, A;/r| w 1.02 x lO'^GeV''. (54) 

As far as the critical value of the Higgs boson mass is concerned, we know explicit expressions for 
the entropy density and the sound speed ( ^5|) in this case. Substitution of the sound speed to the 
equations for the critical point gives the result 

2 \/33-3 2 7-\/33 3rg 

Then the equation for the entropy flux conservation takes the form 



r^a /3 



r^u^iC^ -^,) = -l^-{63 + nVsS) = r^^e " 1 + rgKC^ - 0- (56) 

This equation can be reduced to the fourth order polynomial equation on 1/r. Only one root of this 
equation has physical meaning and satisfies the boundary condition (|48|). The explicit expression 
for r(^^) is rather huge and we do not write it here, but it is remarkable that the accretion problem 
admits of an analytical solution. 

We see from the plots presented in Fig. 13 that the singularities arising in the effective potential 
at ^^ = are not actually realized. In particular, when A2 = 0, the minimal value of the reciprocal 
temperature squared, which it takes on the horizon r = Vg, becomes 

C^|^=^^ = ^ (^1402 + 66\/33 - -\/378 + 66\/33 j « 0.46. (57) 

As for the recommended Higgs boson mass, the numerical analysis shows that ^^ ~ 0.62 at r = r^. 
This corresponds to the relative change of the vacuum expectation value 6r]/r]Q ~ 9.2 x 10^'^. 
Taken on the horizon, the Tolman reciprocal temperature squared ^^ decreases, when the Higgs 
boson mass rriHo increases till mno ~ 232 GeV. The reciprocal temperature squared is equal 
approximately to 0.56 at this value of the Higgs boson mass. Then it increases up to unity at ffiHo- 
After the instability interval, the reciprocal temperature squared increases with mno starting from 
(57). At very large Higgs boson masses, applicability of the perturbation theory and the one-loop 
approximation becomes questionable and so we exclude this region from our consideration. 
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FIG. 4. The dependence of ^'^ on the Schwarzschild coordinate r. Two branches of the solution to the accretion 
equations (Uq) are presented. One branch corresponds to the "infaU" solution, another one being the "outflow" 
solution. These branches intersect at the critical point. The insets depict the moduli of the accretion velocities. 
Only that branch of the general solution which does not vanish on the horizon represents the solution satisfying the 
boundary condition (ffq). This is the infall solution. The case of the recommended Higgs boson mass is presented on 
the left panel and the case of the critical Higgs boson mass m^g is given on the right panel. These plots together 
with the plots in Fig. H provide a complete description of the gravitational mass-shift effect for a Schwarzschild black 
hole. 



V. DISCUSSION 

In this paper we obtained the one-loop effective potential of the Higgs field both for stationary 
and non-stationary gravitational backgrounds. The crucial point was the use of the energy cutoff 
regularization or the normal ordering in prescribing a rigorous meaning to the Hamiltonian of 
the standard model and to the other composite operators. Having imposed physically reasonable 
normalization conditions, we completely specified the effective potential and described thereby the 
vacuum state at the one-loop level. This, in turn, allowed us to obtain a complete description of 
the gravitational mass-shift effect on the Schwarzschild background in a certain approximation. 

In particular, it appeared that the gravitational mass-shift effect is greater for static stars than 
for Schwarzschild black holes. For the non-rotating neutron star with the radius of 2rg this effect 
gives the relative mass-shift 6r]/r]Q « 2.0 x 10^^ (for the recommended value of the Higgs boson 
mass) on the surface of the star, while for the black hole it leads to 6r]/r]Q ~ 9.2 x 10^^ on the 
horizon. A proper generalization of the approach to the non-stationary case allowed us to get rid 
of the divergences of the effective potential on the black hole horizon. 

Besides, the properties of the vacuum proved to be similar to the properties of an ultrarelativistic 
fluid. It possesses the entropy and enthalpy densities, the pressure etc. The entropy density and the 
pressure of the vacuum turn out to be negative, when the Higgs boson mass is less than 263.6 GeV, 
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and they become positive for the Higgs boson masses greater than the critical value 278.2 GeV. 
This implies the existence of a small screening of a gravitating object in the former case and a small 
anti-screening in the latter case. Although these effects are rather small, they become relevant on 
cosmic scales. The vacuum energy density tends to zero as —a'{l)rg/r and the pressure behaves 
like p"{l)r^/2r'^ at sufficiently large distance r from the object. Therefore the energy of the vacuum 
state diverges in the limit of infinite space. Of course, this contribution to the vacuum energy is only 
relevant on the distances less than |cr'(l)|rg/A as the contribution from the cosmological constant 
A dominates above this scale. 

The results of this paper can be generalized in several directions. It would be interesting to 
investigate the properties of the vacuum state for the Kerr-Newman background. A naive substi- 
tution of the Kerr-Newman metric to the effective potential (|23| ) gives rise to the divergence on 



the ergosphere. This problem seems to be resolved by solving the self-consistency condition (37) as 
we have done for the Schwarzschild black hole. Another evident generalization is to include back- 
reaction and derivative corrections to the one-loop effective action, although, to all appearance, 
these corrections are small for macroscopic gravitating objects. It is also interesting to investigate 
the loop corrections to other observables, such as the electron form factors on a curved background 
using the energy cutoff. An evident generalization of the flat spacetime results suggested by the 



form of the effective potential (23) consists in replacement of the massive parameter fj, of the dimen- 
sional regularization (or the cutoff parameter) by its blue-shifted counterpart /u(^^)^^'^ provided 
the derivatives of ^^ are negligible. This will lead to small variations of the coupling constants with 
gravity (see for recent tests, e.g., ||57|| ). However, this guess needs a further exploration. When this 
problem will be solved, the higher loop corrections to the effective potential of the Higgs field on a 
curved background can be studied making use of the energy cutoff regularization. 
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